Nonexponential Relaxations in a Two-Dimensional Electron System in Silicon 
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The relaxations of conductivity have been studied in a strongly disordered two-dimensional (2D) 
electron system in Si after excitation far from equilibrium by a rapid change of carrier density ris at 
low temperatures T. The dramatic and precise dependence of the relaxations on ria and T strongly 
suggests (a) the transition to a glassy phase as T 0, and (b) the Coulomb interactions between 
2D electrons play a dominant role in the observed out-of-equilibrium dynamics. 
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The interplay of strong electronic correlations and dis- 
order is believed to be responsible for a plethora of new 
phenomena occurring in many complex materials in the 
region of the metal-insulator transition (MIT) [1]. Even 
though glassy behavior of electrons is emerging as one of 
the key concepts, electron or Coulomb glasses [2] in gen- 
eral remain largely unexplored. Slow, nonexponential re- 
laxations characteristic of glassy dynamics have been ob- 
served in doped semiconductors [3], strongly disordered 
InOa: films [4-6] , and various granular metals [7] , but the 
regime near the MIT has been even less studied. In many 
complex systems, the MIT is accompanied by changes in 
magnetic or structural symmetry, which complicates the 
situation even further. On the other hand, low-density 
two-dimensional (2D) electron and hole systems in semi- 
conductor heterostructures, where the MIT has been a 
subject of great interest and debate [8], represent partic- 
ularly appealing model systems for studying the effects of 
interactions and disorder in a controlled and systematic 
way. 

Our recent resistance noise measurements on a 2D elec- 
tron system (2DES) in Si [9-11] have shown signatures 
of glassy behavior at densities ris lower than some well- 
defined density Ug, such that rig > Uc {ric - the critical 
density for the MIT). In particular, by reducing n^, it 
was observed that the dynamics suddenly and dramati- 
cally slowed down at Ug , and there was an abrupt change 
to a correlated statistics, consistent with the hierarchical 
picture of glassy dynamics. These features were shown to 
persist even when the 2DES is fully spin polarized [11], 
indicating that the charge degrees of freedom are respon- 
sible for the anomalous noise behavior. Similar results 
were obtained in both highly disordered samples [9] and 
those with relatively low disorder [10, 11] , with rig ^ ric in 
the latter case. In addition, for ris < rig, the conductivity 
(T of the 2DES was found to depend on the cooling pro- 
cedure [9], another characteristic of glassy systems. We 
note that, even though rig is in the metallic regime (i. e. 
(t(T ^ 0) 7^ 0), a is so small that kpl < 1 {kp - Fermi 
wave vector, I - mean free path) [12], which violates the 
Mott limit for the metallic transport in 2D. Hence the 
noise studies [9-11] as well as this work probe the regime 
of strong disorder that encompasses both the insulating 



phase and the unconventional conducting regime where 
kpl < 1. Such "bad metals" include a variety of strongly 
correlated materials with unusual properties [13]. 

Here we report a systematic study of the relaxations 
(i.e. time dependence) of a in the same highly disor- 
dered 2DES in Si after applying a large perturbation at 
different and T. The perturbation consists of a large, 
rapid change of ris, controlled by the gate voltage Vg. For 
ris < t^g, we find that, on short enough time scales, a[t) 
obeys a nonexponential, Ogielski form [14], found in some 
other systems [15] just above the glass transition. The 
associated relaxation time tiow diverges as T — > 0, while 
its density dependence provides strong evidence that the 
observed out-of-equilibrium behavior is dominated by the 
Coulomb interactions between 2D electrons. After a suf- 
ficiently long time t ^ Thigh, c relaxes exponentially to 
its equilibrium value. The corresponding characteristic 
time Thigh ^ oo as T — > 0, suggesting that the glass 
transition temperature Tg = 0. The most peculiar fea- 
ture of the relaxation is that the system equilibrates only 
after it first goes farther away from equilibrium. 

Measurements were carried out on a 2DES in Si metal- 
oxide-semiconductor field-effect transistors (MOSFETs) 
with a large amount of disorder (the peak mobility 
K, 0.06 m^/Vs at 4.2 K with the applied back-gate bias 
of —2 V). The sample dimensions LxW {L- length, W- 
width) were 1 x 90 /im^ and 2 x 50 fiw? . Both samples 
exhibited the same behavior and, in general, they were 
almost identical to the samples used previously in noise 
measurements [9]. The data obtained on a 2 x 50 ^.tc? 
sample are presented in detail below. For this sample, 
n,.(10"cm-2) = 4.31(K,[V] - 6.3); ng(10"cm-2) « 7.5 
and nc(10^^cm^'^) « 4.5, where rig was determined from 
noise and from a{ns, T) measurements on both metal- 
lic and insulating sides (see Refs. 9, 11). The samples and 
the standard ac lock-in technique (typically 13 Hz) that 
was used to measure a have been described in more detail 
elsewhere [9]. 

The experiment consists of the following procedure. 
The sample is cooled from 10 K to the measurement 
temperature T with an initial gate voltage Vg. Then, 
at t = 0, the gate voltage is switched rapidly (within 1 s) 
to a final value , and (j{t, , T) measured. In general. 
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FIG. 1: (Color online) (a) a{t) for = 11 \ 

K(10"cm-2) = 20.26], V/ = 7.4 V ^(lO^cm-^) = 4.74] 
and T = 3.3 K. (b) Experimental protocol: Vg{t) and T{t). 

the results did not depend on any of the following: initial 
temperature, as long as it was > 10 K; Vg] the cooling 
time, which was varied between 30 minutes and 10 hours; 
the time the sample was kept at 10 K (from 5 minutes to 
8 hours); the time the sample spent at the measurement 
T before Vg was changed (from 5 minutes to 8 hours). 
Figure 1 shows a typical experimental run with a{t > 0) 
exhibiting a rapid (< 10 s) initial drop followed by a 
slower relaxation [16]. After going through a minimum, 
a{t) increases and approaches a value (To{VJ' , T). A sub- 
sequent warm-up to 10 K and a cool down to the same 
measurement T = 3.3 K, while keeping the gate voltage 
fixed at , shows that (Jo{V^ ,T) represents the equi- 
librium conductivity corresponding to the given and 
T . It is interesting that, even though initially it drops to 
a value close to ctq, <t first goes away from equilibrium 
before it starts approaching ctq again. 

At the end of the run, the sample is warmed up to 10 K, 
gate voltage changed back to the same , and the exper- 
iment is repeated at a different T for the same . Fig- 
ure 2 shows the relaxations of conductivity, a{t,Vj ,T), 
normalized to the corresponding aoiVj" ,T) for different 
T. The minimum in a shifts to longer times with de- 
creasing T until, at sufficiently low T, it falls out of the 
time window of the measurements. The relaxations ob- 
viously become slower with decreasing T but the shape 
of the relaxation curves remains similar, suggesting that 
they might be collapsed onto a single scaling function. 

First we focus on times just before the minimum in 
a{t). Figure 3(a) shows that the data can be described 
by a scaling form a{t,T)/aoiT) = a{T)git/Tio^{T)), 
such that a(T) oc (rfo^j)^" (a ~ 0.07) and git/now) ^ 
eyi\)[—{t/TiowY], P = 0.3. We note that this stretched 
exponential dependence spans more than 4 orders of 
magnitude in t/riow^ and 2 orders in cr/tTo- At shorter 
times, most easily observable at low T within our ex- 
perimental time window, the stretched exponential func- 
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FIG. 2: (Color online) cr(t) for = 11 V, 1// = 7.4 V, and 
several T, as shown. For 3.2 < r(K) < 4.2, T was changed in 
steps of 0.1 K but some traces have been omitted for clarity. 
Similarly, the data obtained at 0.24 < T{K) < 1 are shown 
in Fig. 3(b). 



tion crosses over to a slower, power-law dependence 
a{t,T)/ao{T) cx [Fig. 3(b)] [17]. In this regime, 
of course, it is also possible to collapse the data but, 
obviously, not in an unambiguous way. Both power- 
law and stretched exponential relaxations are consid- 
ered to be typical signatures of glassy behavior, and 
refiect the existence of a broad distribution of relax- 
ation times. In spin glasses, for example, the so- 
called Ogielski function [14], g{t) cx i"" exp[-(i/r(T))'5], 
describes the relaxations over a wide range of t and 
T above Tg [15]. Our results [Figs. 3(a), (b)] sug- 
gest that the Ogielski function might be a good can- 
didate to describe also our data over the entire t and 
T range. Indeed, Fig. 3(c) shows that all the data are 
consistent with the Ogielski relaxation a{t,T)/aa{T) cx 
{Tiow)~"'{t/Tion])~"' exp[-{t/Tiow)~^] ovcr about 25 or- 
ders of magnitude in t/riow Since it was impractical 
to perform measurements below 1 K long enough to ob- 
tain a substantial overlap between different T curves, we 
limit our analysis of Tiom (T) to the T > 1 K range, where 
the data collapse well. Even within this relatively nar- 
row range of T, the scaling parameter tiow varies over 7 
orders of magnitude [Fig. 3(d)], and the Arrhenius func- 
tion, 1/t,ou, = A:oexp(-£:a/r) (fco 10 s-\ Ea ~ 19 K 
for Vg^ = 7.4 V), provides an excellent fit to the data [18]. 
The experiment was repeated for different values of 
, allowing us to map out the density dependence of 
the exponents a and (3, and of the parameter tiow The 
results were obtained either by fitting the data to the 
appropriate functional form at a fixed T [as in Figs. 3(a) 
and (b)] or from the scaling of the different T data and 
the fit to the Ogielski function [as in Fig. 3(c)]. Figure 4 
shows that the two methods yield very similar results. 
In particular, we find that a decreases with n^, and van- 
ishes at ns(10^^cm~^) « 7.5 « Ug, i. e. exactly where 
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FIG. 3: (Color online) (a) Scaling of the Fig. 2 data for 
t just below the minimum in u{t). The data obey a 
stretched exponential dependence on t; the dotted line is 
a fit. The prefactor a{T) cx (tio^)'", a 0.07. (b) 
At short t, relaxations follow a power law a{t,T)/ao{T) oc 
t'" [17]. The dashed lines are fits with slopes a — 0.07 ± 
0.01. (c) While the lowest T data clearly deviate from the 
stretched exponential dependence (dashed line), all the data 
are consistent with the Ogielski relaxation, a{t,T)/ao{T) oc 
iriow)~'^it/Tiow)~" exp[-{t/Tio-w)~'^] (dotted line). The data 
have been collapsed with respect to the 2.4 K curve, (d) Scal- 
ing parameter r;o„ as a function of T. The dashed line is 
a fit with a slope equal to an activation energy Ea ~ 19 K; 
r,„„(2.4 K) ~ 500 s. 



the onset of glassy dynamics was found to take place from 
the noise measurements [9]. At the same time, the ex- 
ponent (3 grows with increasing Hs, indicating that the 
relaxations become faster, as expected. For all Ug in this 
range, 1/tiow = ko{ns) exp{-Ea/T) with Ea ^ 20 K 
independent of rig, but the prefactor has a very strong 

1 /2 

density dependence, ko{ns) oc exp(— 7ns ) [Fig. 4(b)]. 

Since the ratio of the Fermi energy to Coulomb energy 

1/2 

in 2D systems l/rg = Ep/U oc ng [8, 21], our result 
strongly suggests that Coulomb interactions between 2D 
electrons play a dominant role in the observed slow dy- 
namics. As rig is increased further, beyond rig, there are 
still some visible relaxations, albeit with an amplitude 
(t/(To that is too small to make reliable fits. In fact, the 
amplitude of the relaxations decreases with Ug , and van- 
ishes when ns(10^^cm~^) 30, where ctq ~ 3e^/h, i. e. 
kpl ~ 1. Therefore, all the phenomena discussed in this 
paper take place in the regime of strong disorder. We 
recall that, in general, electron-electron interactions are 
enhanced in the presence of disorder [22, 23]. 

Finally, at times above the minimum in a{t) (Fig. 2), 
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FIG. 4: (Color online) (a) The exponents a and [3 vs. Us- 
Dashed lines guide the eye. (b) logrjom plotted vs. n'J'^ . The 
values of riam correspond to T = 3 K. The dashed line is a fit 
with a slope (log 6)7 4.5(10~"cm^)^''^ 




FIG. 5: (Color online) (a) Scaling of the Fig. 2 data at long 
t, i. e. above the minimum in (j{t). The data have been col- 
lapsed with respect to the 5 K curve, (b) The scaling function, 
which describes the approach to equilibrium value ao at long 
times. The dashed lines are fits with slopes k\ ~ 2.89 x 10~^ 
and ^2 ~ 4.37 X 10~^ for the faster and slower exponential pro- 
cesses, respectively. (The corresponding relaxation times are 
Thigh/ki and Thigh/k2.) Inset: T-dependence of Thigh/ki for 
different Vg and V/ as shown on the plot. The dashed line is 
a fit with the slope equal to an activation energy Ea ~ 57 K. 

relaxations at different T approach the corresponding 
equilibrium values cro(T). Here all the data can be de- 
scribed by a scaling form a{t,T)/ao{T) = f {t / Thigh{T)) 
[Fig. 5(a)], where f{t/ Thigh) is a simple exponential func- 
tion with the characteristic relaxation time Thigh /ki. For 
some Us [see Fig. 5(b)], another slower, simple expo- 
nential process also seems to exist with the same T- 
dependence of its characteristic relaxation time Thigh /k2, 
but its presence does not seem to depend on in any 
systematic way. We find that, unlike ti^w, Thigh/ ki 
(and Thigh /k2, if present) does not depend on iig in 
the regime studied. On the other hand, ki/Tfiig}i — 
(2.7x 10^ s-i) exp(-£;A/r) with Ea^57K independent 
of Us (Fig. 5(b) inset). 

Several striking features of our data stand out. 

(1) At long times, the relaxation is exponential with 
a characteristic time Thigh /ki that diverges exponen- 
tially as T — > 0. This means that, at any finite T, the 
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2DES will reach equilibrium after a sufficiently long time 
t > {Thigh /ki). For example, Tugh/ki ~ 30 s at T = 5 K 
but, already at T = 1 K, Thigh /k\ ^ 10^^ years, which 
exceeds the age of the Universe by several orders of mag- 
nitude. Therefore, even though, strictly speaking, the 
system appears to be glassy only at T = 0, at low enough 
T the dynamics is glassy on all experimentally accessible 
time scales. It is interesting that the recent study of the 
2D Coulomb glass model has also found [24] an expo- 
nential divergence of the equilibration time as T — »■ 0, 
signaling a transition at Tg = 0. 

(2) In the glassy regime, i. e. on short enough time 
scales t < {viigh/ki), the relaxations are strongly nonex- 
ponential. They obey the Ogielski form a{t,T)/ao{T) oc 
t-"exp[-(t/rio^)-'5] with < a < 0.4 and 0.2 < /3 < 

1 /2 

0.45, where tiow oc exp(7ns ) cxp{Ea/T), Ea ~ 20 K. 
Since tiow diverges as T — > 0, the relaxations attain a 
pure power law form oc at T = 0. We note that the 
Ogielski form of the relaxation and the divergence of tiow 
are consistent with the general scaling arguments [25] 
near a continuous phase transition occurring at Tg = 0. 
Similar scaling is observed in spin glasses above Tg [15]. 

The very pronounced and precise dependence of Tiow 
on rig strongly suggests that Coulomb interactions be- 
tween 2D electrons play a dominant role in the observed 
out-of-equilibrium dynamics. In InO^; films, the most ex- 
tensively studied example of an electron glass [4-6], the 
dependence of the relaxation time on the carrier density 
has been interpreted as evidence that the slow dynamics 
is due to electronic relaxation rather than glassiness of 
the extrinsic degrees of freedom [5] . 

(3) The 2DES equilibrates only after it first goes far- 
ther away from equilibrium. This overshooting of a man- 
ifests itself as cither a minimum [Fig. 1(a)] or a maximum 
in a{t) depending on the direction of Vg change. For ex- 
ample, ifVg < Vj , (T{t) first increases, overshoots uo and 
reaches a maximum, and then decreases exponentially 
to (To. An analogous overshooting during relaxation is 
known to occur in orientational [26] and spin glasses [27] 
but has not been observed in InOj; films [6] . We note an 
important difference between our experiment and those 
on InOx films. In InO^, the change in the carrier density 
due to variations in Vg is typically 1%, and the system 
remains deep in the insulating state. In our case, Us is 
changed up to a factor of 7, and the 2DES may go from 
the conducting to the insulating regime. Recent numeri- 
cal work has shown [28] that such a strong perturbation 
can lead to nontrivial dynamics and, in particular, to a 
"roundabout" relaxation, where system equilibrates only 
after it once goes farther away from equilibrium. 

In summary, we have studied the relaxations of tr in a 
strongly disordered 2DES after excitation far from equi- 
librium by a rapid change of at low T. The data 
strongly suggest (a) the transition to a glassy phase as 
T ^ 0, and (b) the Coulomb interactions between 2D 
electrons are primarily responsible for the observed out- 



of-equilibrium dynamics. Further work is needed to iden- 
tify the microscopic details of the relaxation phenomena. 
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